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Abstract 

The form factors relevant to B q —* D*(J P = l~)£v (q = s,d,u) 
decays are calculated in the framework of the three point QCD sum 
rules approach. The heavy quark effective theory prediction of the 
form factors as well as 1/mb corrections to those form factors are ob- 
tained. A comparison of the results for the ratio of form factors at 
zero recoil limit and other values of q 2 with the predictions of the sub- 
leading Isgur-Wise form factor application for B — ► D*lv is presented. 
The total decay width and branching ratio for these decays are also 
evaluated using the q 2 dependencies of these form factors. The results 
are in good agreement with the constituent quark meson model and 
existing experimental data. The q = s case can also be detected at 
LHC in the near future. 
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1 Introduction 



Semileptonic pseudoscalar B q decays are crucial tools to restrict the Stan- 
dard Model (SM) parameters and search for new physics beyond the SM. 
These decays provide possibility to calculate the elements of the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix, leptonic decay constants as well as the 
origin of the CP violation. 

When LHC begins to operate, a large number of B q mesons will be pro- 
duced. This will provide experimental framework to check the B q decay 
channels. An important class of B s and entire B u ^ decays occur via the 
b quark decays. Among the b decays, b — > c transition plays a significant 
role, because this transition is the most dominant transition among b decays. 
Some of the B decay channels could be B q — ► D*lv (q = s,d,u) via b —>■ c 
transition. These decays could give useful information about the structure 
of the vector D* mesons. The observation of two narrow resonances with 
charm and strangeness, _D s j(2317) in the D s tc° invariant mass distribution 
[I]-©, and L> sJ (2460) in the D*ir° and D sl mass distribution [21 EJ m E3 E] , 
has raised discussions about the nature of these states and their quark con- 
tents EUj. Analysis of the £> So (2317) -> D*~f, D sJ (2460) -> D* sl and 
D s j(2A60) — > D SQ (2317)7 indicates that the quark content of these mesons 
is probably cs [TT] . 

Form factors are central objects in studying of the semileptonic B q — > 
D*lv decays. For the calculation of these form factors, we need reliable non- 
perturbative approaches. Among all non-perturbative models, the QCD sum 
rules has received especial attention since this model is based on the QCD 
Lagrangian. QCD sum rules is a framework which connects hadronic pa- 
rameters with QCD parameters. In this method, hadrons are represented 
by their interpolating currents taken at large virtualities. The correlation 
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function is calculated in hadrons and quark-gluon languages. The physical 
quantities are determined by matching these two representations of correla- 
tors . The application of sum rules has been extended remarkably during the 
past twenty years and applied for wide variety of problems ( For review see 
for example [12J). 

The aim of this paper is to analyze the semileptonic B q — > D*lu decays 
using three point QCD sum rules method. Note that, this problem has been 
studied for B q — > D*lv (q = s, d, u) in constituent quark meson (CQM) model 
in [13] and for q = d 1 u{B Q 1 B ± ) cases in experiment [14] . The application 
of subleading Isgur-Wise form factor for B — > D*lv at heavy quark effective 
theory (HQET) is calculated in [15] (see also [16j[T7]). Present work takes 
into account the SU(3) symmetry breaking and could be considered as an 
extension of the form factors of D — > K*ev presented in [18] . 

The paper is organized as fallow: In section II, sum rules expressions 
for form factors relevant to these decays and their HQET limit and l/m^ 
corrections are obtained. The numerical analysis for form factors and their 
HQET limit at zero recoil and other values of y, conclusion and comparison 
of our results with the other approaches are presented in section III. 

2 Sum rules for the B q —> D* q iu transition form 
factors 

The B q — > D* transitions occur via the b — > c transition at the quark level. 
At this level, the matrix element for this transition is given by: 

G 

M q = -^V cb V Tm (1 - l5 )l c 7m (1 - T5 )6. (1) 

To derive the matrix elements for B q — > D*lv decays, it is necessary to 
sandwich Eq. ([1]) between initial and final meson states. The amplitude of 
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the B q — > -D*/f decays can be written as follows: 

M = ^V cb V 7^(1 - 75 )/ < £>*(?/,£) | c 7m (1 - 75 )& I S,(p) > • (2) 

The aim is to calculate the matrix element < D*(p',e) | C7 M (1 — 75)6 | 
B q (p) > appearing in Eq. ([2]). Both the vector and the axial vector part of 
c 7^(1 — 75)6 contribute to the matrix element stated above. Considering 
Lorentz and parity invariances, this matrix element can be parameterized in 
terms of the form factors below: 

< D* q (p', e) I c 7 ,6 I B 9 (p) >= i fv{ f ] ^ ^Y, (3) 

(m Bq +m D *) 



< D * q (p', £ ) I C7 M 75& I B q (p) > = i [f (q 2 )(m Bq + m D *)e 
U{q2) (^)P„ - , Uq2) & p)*> 



(4) 



(m Bq + m D *) (m Bq + m D 

where fvio 2 )-, foiQ 2 )-, /+(? 2 ) an d f-{<f) are the transition form factors and 
Pfi = (p + p')m> % = (P ~ p')v m or der to calculate these form factors, 
the QCD sum rules method is applied. Initially the following correlator is 
considered: 

nJ2V,P'W) = J d^ye-^'y < I T[J„ Dq (y)jV-> A (0)J Bq (x)) | >, 

(5) 

where J U D*(y) = QluC and J Bq (x) = try^q are the interpolating currents of D* 
and B q mesons, respectively and = c^^b and J A = 07^756 are vector and 
axial vector transition currents . Two complete sets of intermediate states 
with the same quantum numbers as the currents Jo* and J Bq are inserted to 
calculate the phenomenological part of the correlation function given in Eq. 
(JSj). After standard calculations, the following equation is obtained: 



< I r D , | D* q (p',e) >< D* q (p',e) \ J^ A | B q (p) >< B q (p) I J Bq 1 > 

(p /2 -m 2 D ,){j> 2 -m% g ) 

(6) 

where • ■ • represents contributions coming from higher states and continuum. 
The matrix elements in Eq. ([6]) are defined as: 

fB m \ 

< | r D « | D* q (p', e) >= f D * q m D *e» , < B q (p) \ J B J >= -i q B \ 

(7) 

where Jd* and f Bq are the leptonic decay constants of D* and B q mesons, 
respectively. Using Eq. ([3]), Eq. (TjJ and Eq. ([7]) and performing summation 
over the polarization of the D* meson in Eq. (jH]) the equation below are 
derived: 



(m b + m q ) {p' 2 - m 2 D *)(p 2 - m 2 Bq ) 

(m Bq + m D *) (m Bq + m D * 
+ excited states, 



x 



[-fog^u(m Bq + m D *) + + ^ - + 



jtV ( 2 12 2x a tfi fB g m B q fp* q m D' q 



(m 6 + m s )(m Bq + m D *) (p' 2 - m 2 D ,)(p 2 - m% a 
+ excited states. 

From the QCD (theoretical) sides, Ii^ y {p 2 ,p 12 , q 2 ) can also be calculated 
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qj 



by the help of OPE in the deep space-like region where p 2 <C (wi(, + m, 
and p' 2 <C (m c + m q ) 2 . The theoretical part of the correlation function is 
calculated by means of OPE, and up to operators having dimension d = 6, 
it is determined by the bare-loop (Fig. 1 a) and the power corrections (Fig. 
1 b, c, d) from the operators with d = 3, < ipip >, d = 4, m s < iptp >, 
d = 5, m\ < ipip > and d — 6, < ipip^ip >■ The d = 6 operator is 
ignored in the calculations. To calculate the bare-loop contribution, the 
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Figure 1: Feynman diagrams for B q — > D*lv (q = s,d,u) transitions. 

double dispersion representation for the coefficients of corresponding Lorentz 
structures appearing in the correlation function are used: 

jjper _ _ _ 1— f ds' ( ds- — ^4rr7 \ 9 - )0 x + subtraction terms. (9) 
(2tt) 2 J J (s - p 2 )(s' - p' 2 ) 

The spectral densities p«(s, s', q 2 ) can be calculated from the usual Feynman 
integral with the help of Cutkosky rules, i.e., by replacing the quark prop- 
agators with Dirac delta functions: p 2^ m 2 — > —2tt5(p 2 — m 2 ), which implies 
that all quarks are real. After long and straightforward calculations for the 
corresponding spectral densities the following expressions are obtained: 



p v (s,s',q) = 4N c I (s,s',q)[(m b -m q )A+(m c -m q )B - m q ] , 



p (s,s',q 2 ) = -2N c I (s,s',q 2 ' 



2rr? — 2m 2 {m c + m h ) 



+ m q (q 2 + s + s' - 2m b m c ) + [q 2 (m b - m q 
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+ s(3m q — 2m c — m b ) + s'{m q — m h )]A + [q 2 {m c — m q ) 
+ s(m q — m L ) + s'{3m q — 2m b — m c )}B + A(m b — m s )C 



p + (s,s',q 2 ) = 2NJ (s,s',q 2 '' 



m q + {3m q — nib) A + (jn q — m c )B 



+ 2(m g + m b )D + 2(m q -m b )E 



p-(s,s,q 2 ) = 2NJ (s,s',q 2 ) 



- m q + (m q + m b )A - (m q + m c )B 



+ 2(m q - m b )D + 2(m b - m q )E 



where 
Io(s, s',q 2 ) 



1 



4A 1 /2( S? s / >g 2)' 

A(a, fo,c) = a 2 + 6 2 + c 2 - 2ac - 26c - 2ab, 

1 



A = 



(s' + s - q 2 ) 2 - Ass' 



2m 2 b + q +s- s)s' 



B 



+ m 2 (q 2 - s + s') + ml(-q 2 + s + s' 



m 2 (q 2 + s - s') 



C = 



(s' + s - q 2 ) 2 - Ass' 
+ (-2m 2 c + q 2 - s + s')s + m 2 b {-q 2 + s + s) 
1 



4 , 4 / 

m c s + m b s 



2[(s' + s- q 2 ) 2 - Ass'] 
+ Q 2 [ m t + m g(<? 2 — s — s') + ss 1 + m l m c( ( l 2 — s — s r ) 

- (q + s - s ') s ' - m g(<? 2 — s + s') 

— mlm 2 (q + s — s') + s(g — s + s') 



(10) 



6 



D 



E 



mt[q 4 - 2q 2 (s - 2s') 



[(s' + s- q 2 ) 2 - Ass'] 2 
+ [6m 4 b + q + q\As - 2s 1 ) + (s - s') 2 - 6m 2 b (q 2 + s - s')]s' 2 

+ m\[q + s 2 + 4ss' + s' 2 -2q\s + a')] 

- 2m\s\-2q + (s - s') 2 + 3m 2 b (q 2 - s + s') + q 2 (s + s')] 

- 2m 2 c m 2 q (q 2 + s 2 + ss' - 2s' 2 + q \-2s + s')) 

+ s '[q 4 + q 2 s — 2s 2 — 2q 2 s + ss' + s' 2 + 3m 2 (— q 2 + s + s')\ 



o 4 4 , 2 6 42 24 42 

2m q + m q — m q s — m q s — m s 



.2, 4 



[s' + s - q 2 ) 2 - 4ss'] 2 

o2o o q A z / 2 4 / r\ 4- / 

— m ? g s + m g s — m g g s — m^g s + 2m q ss 

2 ^* / ^ / 2 2/ ^ 2 / / 

+ 6m g ss + 2g ss — m s s — g s s — s s 



i o 4/ 4/2 2 

+ 3m 6 (g — s + s Js — m s — m g s 



./2 



2 /2 

m q ss 



- q ss' 2 + 2s 2 s' 2 + m 2 q s' 3 - ss' 3 - 3m^s(-g + s + s') 

- 2m 2 c m 2 q [q 4 - 2s 2 + q 2 (s - 2s') + ss' + s' 2 )] 
+ s[q 2 + s 2 + ss' - 2s' 2 + q 2 (-2s + s')] 

+ 2ml{— m 2 (q 4 — 2q 2 s + s 2 + q ' s + ss — 2s' 2 ) 

- s'(q 4 + q 2 s -2s 2 - q 2 s' + ss' + s' 2 ) 

+ m^[g 4 + s 2 + 4ss' + s /2 -2g 2 (s + s / )]} 



(11) 



The subscripts V, and ± correspond to the coefficients of the structures 
proportional to ie^ a (3p' a p^ SV and \{p^Pv±p'i 1 Pv)- l respectively. In Eq. (JTUJ) 
N c = 3 is the number of colors. 

The integration region for the perturbative contribution in Eq. is 
determined from the condition that arguments of the three S functions must 
vanish simultaneously. The physical region in s and s' plane is described by 
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the following inequalities: 



2ss' + (s + s'- g 2 )(m 2 b - s - m 2 ) + (m| - m 2 )2s 

A 1 /2( m 2 )S)m 2) A l/2( S!S / ;g 2) -+• \ ) 

From this inequalities, we calculate s in terms of s' in order to put to the 
lower limit of integration over s. For the contribution of power corrections, 
i.e., the contributions of operators with dimensions d — 3, 4 and 5, the 
following results were derived: 



+ + = - 2 <- qq > 



1 

1 



m 2 (m 2 -2m 2 ) 



r[-3m (m b +m c -q ) 



+ 

/ (3) + /o (4) + / (5) = \<m> 



+ ml(ml + m b m c + ml - q 2 )} 

1 1 

" - ^v m '( m o - 2m 2 ,) 

1 2 

^v( 2m o - 3m 6 m,) + — 



x (m 2 + 2m b m c + m 2 c — q 2 ) 

'm\ + 2m b m c + m 2 — g 2 ) 



3r 2 r' 2 

\m\(m\ + 3m b m c - q 2 ) 



x [-3m 2 (m 2 + m 2 - q 2 ) + m\{m 2 h + m b m c + m 2 - g 2 )] 
1 



3rr' 2 

+ 3(m c — m q )m q (m 2 + 2mjm c + m 2 — g 2 )] 
- -^jm 2 (mo - 2m 2 ) (m 2 + 2m b m c + m 2 - g 2 ) 

+ 7^y[-3(m 6 - m q )m q (m 2 b + 2m b m c + m 2 c - q 2 ) 

+ ml(m c + 3m b m c — q 2 )} 

+ - — :(4mQ + 6m 2 + \2m b m c + 6m 2 
3rr 
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— 3m b m q + 3m c m q — 6m — 6q 



if + /r + /- 



(4) 



f(5) 



J < 9? > 



2mi 



+ 



1 



[-3m^(mj + m 2 — 



+ m (m b + m b m c + m c 

m c m q <qq> + -^-m 2 b (ml - 2m 2 q ) 



+ 
+ 



1 



; (3) + f ,y + 



(4) 



e(5) 



3r 2 r' 
1 



-4m 2 , + 3m q (m b + 2m ? )] 



3rr' 



< qq > 



2/ 2 



2m„ 



^V2"[- 3m ?( m ' + m 
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i 2/ 2 , , 2 2\i 

+ m Q {m b + m b m c + m c - g )J 

1 1 

- — -m c m q 
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m 2 b (m 2 -2m 2 q ) 



+ 



1 2 
—m q {-m b + 2m g ) H 



(13) 

where r = p 2 — m\ and r' = p /2 — m 2 . Here we should mentioned 
that, considering the definition of double dispersion relation in Eq. (Q and 
parametrization of the form factors and the coefficient of selected structures, 
with the changes: 1) b — > c and c — > s, 2) set the m q — > and 3) ignore 



the terms ~ m 2 , the Eqs. fflOl 1131) reduce to the expressions for the spectral 
densities and quark condensate contributions up to 5 mass dimensions for 
the form factors fy, fo and /+ presented in the appendix A of [18] which 
describes the form factors of D —>■ K*ev. 

By equating the phenomenological expression given in Eq. (jSJ) and the 
OPE expression given by Eqs. (11 OH 131) . and applying double Borel transfor- 



J2 



Ml) 



m 



mations with respect to the variables p and p (p — > M 2 , p' 
order to suppress the contributions of higher states and continuum, the QCD 
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sum rules for the form factors fv, fo, f+ and /_ are obtained: 

2 (m b + m q ) r] m 2 Bq /M?+m 2 D jMi 
fB t m% q f D *m D , 

' r " da' r ds Pl (s,s',q 2 )e- s / M ?- s '/ M " 



(2ir) 2 J(m c +m a Y Jf(s') 

+B(fP + ft ] + if)] 



(14) 



where i — V, and ±, and B denotes the double Borel transformation op- 
erator and 7? = 777n + mn. for i = V, ± and 77 = 7 for 7 = are 

' D l u q ' ' m B +m D * 

considered. Here k = +1 for i = ± and k = — 1 for 7 = and V. In Eq. (I14p . 
in order to subtract the contributions of the higher states and the continuum, 
the quark-hadron duality assumption is used, i.e., it is assumed that 

p higherstate S ^ ^ = p OPE ^ _ _ ^ ^ 

In calculations the following rule for the double Borel transformations is used: 

B—— -> (-l) m+n — — e -™ 2 b /Ml e -mi/MZ 1 

r m r' n y ' T{m)T(n) (M 1 2 ) m - 1 (M|)™- 1 ' 1 ' 

Here, we should mention that the contribution of higher dimensions are 
proportional to the powers of the inverse of the heavy quark masses, so this 
contributions are suppressed. 

Next, we present the infinite heavy quark mass limit of the form factors 
for B q — > D*lv transitions. In HQET, the following procedure are used (see 
[TH [201 [21] ) . First, we use the following parametrization: 

777r + 777n* — Q 2 

y = vv= a™. m ~ (17) 
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where v and v' are the four-velocities of the initial and final meson states, 
respectively and y = 1 is so called zero recoil limit. Next, we try to find the 
y dependent expressions of the form factors by taking mj — > oo, m c = ^7=, 
where z is given by ^fz = y + \Jy 2 — 1 and setting the mass of light quarks 
to zero. In this limit the Borel parameters take the form M 2 = 2T\m b and 



Ml 



2T2in c where 7\ and T2 are the new Borel parameters. 



The new continuum thresholds z/ , and u' take the following forms in this 
limit 

~ ? '" s 'o- m l ( lg ) 



s - mi 
= — , 



z/ n = 



m b m c 
and the new integration variables are defined as: 



m b 

The leptonic decay constants are rescaled: 
/s„ = y/rnbf B „, Jd* 



m c f. 



(19) 



(20) 



After the standard calculations, we obtain the y-dependent expressions of 
the form factors as follows: 



, (1 + v/i) 



fo 



MfD-fB^ 1 '* 

— - = / dv \ dv\v + u')e~ w ^~^9(2yuu' - v 2 - u' 2 ) 



i- 1/0 rv. 



+ 16 < qq > 



mU 1 1 1 _L_ 1, 

8 \2T 2 + 2T| + 3TxT 2 { + ^ + z' 



(21) 



— - — e T i T 2 ; < - r / ctz/ / di/(i/ + i/)e 2T i 2T 2 

16/u;/B,(l + >/i) l^v^Tio 7o 



0(2yi/i/ - z/ 2 - z/' 2 ) + 



< qq > yfz 



1 J_ _1_ N 

^2 + 2^ + 
1 



9 , 1 1 ,\ / 1 4 1 Vi N 



(22) 
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A 



+ ^ e ( i+i>< 



96/^/5,^/4 



— ^-^== f° ^ f ° dv\v + z/)e~*-^(2yzV - z/ 2 - z/ 2 ) 



16 < qq > 



1 + 



+ 



+ 



1 + 



2T X 2 2T| 3T\T 2 V ^ 



+ 



2; z 



(23) 



f- 



' 1 + v^) 



96/0-/5^1/4 



9 /'"• /•":, _ _ 

— / dv \ dv\v + v')e ^ ^ 6(2yvv' - v 2 - u' 2 ) 



+ 16 < gg > 



1 _ H!° 



1 1 111^ 
2Tf + 27? + mT 2 ^ + 7l + ? 



(24) 



where A = rriB q — vn^ and A = mp* — m c . 

At the end of this section, we would like to present ^- corrections for 
the form factors in Eqs. (|2TT) -( l24l) using sub leading Isgur-Wise form factors 
similar to [15] (see also [20j[22]). These corrections are given as: 



f(!/ m 6) 
'V 



ft 



(l/m b ) 



m B + m* D f A A 1 



(l/m b ) 



(l/m 6 ) 



1 

2 



rriB + m* D [2mby + l nib 
Jv ' 



y + r 



(l/m b ) 



(25) 



where the explicit expressions for Pi(y) functions are given in [15]. The value 
of those functions at zero recoil limit (y — 1) are given as 



Pl (l) = p 2 (l)=0, p 3 (l)-0, p 4 (l)^i 



(26) 
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3 Numerical analysis 



This section is devoted by numerical analysis for the form factors fv{(f)i 
•/b(<? 2 ), i+(<? 2 ) an d f-(q 2 )- From sum rule expressions of these form factors it 
is clear that the condensates, leptonic decay constants of B q and D* mesons, 
continuum thresholds sq and s' and Borel parameters Mf and Mf are the 
main input parameters. In the numerical analysis the values of the con- 
densates are chosen at a fixed renormalization scale of about 1 GeV. The 
values of the condensates are [23] : < uu >=< dd >= —(240 ±10 MeV) 3 , 
< ss >= (0.8 ± 0.2) < uu > and tuq = 0.8 GeV 2 . The quark masses 
are taken to be m c (u = m c ) = 1.275 ± 0.015 GeV, m s = 95 ± 25 MeV, 
m u = (1.5 - 3) MeV, m d ~ (3 - 5) MeV [H] and m b = (4.7 ± 0.1) GeV 
|23j . The mesons masses are chosen to be mo* = 2.112 GeV , mo* — 
2.007 GeV , m D * = 2.010 GeV, m Bs = 5.3 GeV, m Bd = 5.2794 Gel/and 
m Bu = 5.2790 Gey[14j. For the values of the leptonic decay constants of B q 
and D* mesons the results obtained from two-point QCD analysis are used: 
f Bs = 0.209 ± 38 GeV [12j, f D . = 0.266 ± 0.032 GeV\H\- For the others 
f BM = 0.14 ± 0.01 GeV and f D . = 0.23 ± 0.02 GeVUM are selected. The 
threshold parameters So and s' are also determined from the two-point QCD 
sum rules: s = (35 ±2) GeV 2 [24J and s' = (6-8) GeV 2 [11]. The Borel pa- 
rameters M 2 and Mf are not physical quantities, hence form factors should 
not depend on them. The reliable regions for the Borel parameters M 2 and 
Mf can be determined by requiring that both the continuum contribution 
and the contribution of the operator with the highest dimension be small. 
As a result of the above-mentioned requirements, the working regions are 
determined to be 10 GeV 2 < M 2 < 25 GeV 2 and 4 GeV 2 < Mf < 10 GeV 2 . 

To determine the decay width of B q — > D*lu, the q 2 dependence of the 
form factors /y(g 2 ), fo(q 2 ), /+(<? 2 ) and f-(q 2 ) in the whole physical region 
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mf < q 2 < { m B q — m D*) 2 ar e needed. The value of the form factors at q 2 = 
are given in Table 1. 







B d -> 


# u -> 


Mo) 


0.36 ±0.08 


0.47 ±0.13 


0.46 ±0.13 


/o(0) 


0.17 ±0.03 


0.24 ±0.05 


0.24 ±0.05 


/+(0) 


0.11 ±0.02 


0.14 ±0.025 


0.13 ±0.025 


/-(o) 


-0.13 ±0.03 


-0.16 ±0.04 


-0.15 ±0.04 



Table 1: The value of the form factors at q 2 = 



The q 2 dependence of the form factors can be calculated from QCD sum 
rules (for details, see [Hl[25]). To obtain the q 2 dependent expressions of the 
form factors from QCD sum rules, q 2 should be stay approximately 1 GeV 2 
below the perturbative cut, i.e., up to 10 GeV 2 . Our sum rules, also, are 
truncated at ~ 10 GeV 2 , but in the interval < q 2 < 10 GeV 2 we can 
trust the sum rules. For the reliability of the sum rules in the full physical 
region, the parametrization of the form factors were identified such that in 
the region < q 2 < 10 GeV 2 , these parameterizations coincide with the sum 
rules prediction. Figs. [21 El HI and [5] show the dependence of the form factors 
/V(<? 2 )? foil 2 )) f+(l 2 ) an d f-(l 2 ) on q 2 . To find the extrapolation of the form 
factors, we choose the following two fit functions. 

n 



htf) = rr ^Ifi - 3 ^4 > (27) 

1 + aq ± pq 2 + 7<r + Ag 4 
where q = q 2 /rn 2 B[i . The values of the parameters fi(0), a, j3, 7, and A are 
given in Tables 2, 3 and 4. 

ii) h 

ft{q2) = (l 2 -mh) + (q 2 -m%Y ^ 
14 



The values for a, b and m 2 it are given in Tables 5, 6 and 7. For details about 
the fit parametrization (ii) which is theoretically more reliable and some other 
fit functions see [26l 127] . These two parameterizations coincide well with the 
sum rules predictions in the whole physical region < q 2 < 10 GeV 2 and also 
for q 2 < region. For higher q 2 , starting from the upper limit of the physical 
region the two fit functions deviate from each other and this behavior is 
almost the same for all form factors. As an example, we present the deviation 
of above mentioned fit functions in Fig. [61 From this figure, we see that in 
the outside of the physical region the fit (i) growthes more rapidly than fit 
(ii). The fit parametrization (ii) depicts that the mg« pole exists outside the 
allowed physical region and related to that one could calculate the hadronic 
parameters such as Qbb*d* (see [26| 128]). 





f(0) 


a 


P 


7 


A 


fv 


0.38 


-2.53 


2.77 


-2.41 


0.03 


fa 


0.18 


-1.77 


0.98 


-0.23 


-3.50 


u 


0.12 


-2.90 


3.66 


-3.72 


-1.69 


f- 


-0.15 


-2.63 


2.72 


-0.99 


-6.48 



Table 2: Parameters appearing in the fit function (i) for form factors of the 
B s -> D*(2112)lv at M 2 = 19 GeV 2 , M 2 = 5 GeV 2 . 





f(0) 


a 


P 


7 


A 


fv 


0.46 


-2.90 


2.99 


0.67 


-5.04 


fo 


0.24 


-0.21 


2.19 


-1.68 


-2.15 


u 


0.13 


-4.21 


9.52 


-16.86 


12.97 


f- 


-0.15 


-3.93 


-8.03 


-13.48 


9.15 



Table 3: Parameters appearing in the fit function (i) for form factors of the 
B u -> D* u (2007)£is at M 2 = 19 GeV 2 , M 2 = 5 GeV 2 . 
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f(0) 


a 


P 


7 


A 


fv 


0.47 


-3.08 


4.83 


-5.95 


2.95 


fo 


0.24 


-2.20 


2.18 


-1.83 


-1.90 


u 


0.14 


-4.13 


8.99 


-15.10 


10.65 


f- 


-0.16 


-3.87 


7.73 


-12.71 


8.26 



Table 4: Parameters appearing in the fit function (i) for form factors of the 
B d -> £>2(2010)&/ at M 2 = 19 GeV 2 , M 2 = 5 GeV 2 . 





a 


b 


m % 


fv 


55.03 


-54.30 


23.18 


fo 


1.43 


-4.32 


18.80 


f+ 


1.14 


-2.57 


14.88 


f- 


-2.80 


3.43 


14.60 



Table 5: Parameters appearing in the fit function (ii) for form factors of the 



B s -> D*(2112)£u at M 2 



19 Gel/ 2 , Mf 



5 GeV 2 . 



In deriving the numerical values for the ratio of the form factors at HQET 
limit, we take the value of the A and A obtained from two point sum rules, 
A = 0.62GeV [29J and A = 0.86GeV[3D]. The following relations are defined 
for the ratio of the form factors, 



#1(2) [3] 
#4(5) 
Rr 



1 - 



(f 



{m B + m D *) 2 

,2 



q 



(m B + m D ») 2 

-,2 



1 - 



q 



(m B + m D ») 2 



fv(+)l-](y) 
fo(y) ' 

f+(-)(y) 
fv(y) ' 

f-(y) 
U(vY 



(29) 



The numerical values of the above mentioned ratios and a comparison 
of our results with the predictions of [12] which presents the application of 
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a 


b 


m )u 


fv 


118.69 


-108.48 


23.43 


fo 


4.54 


-5.12 


20.74 


u 


7.79 


-5.84 


14.57 


f- 


-6.72 


5.46 


14.02 



Table 6: Parameters appearing in the fit function (ii) for form factors of the 
B u -> D* u (2007)£is at M\ = 19 Gel/ 2 , Mf = 5 GeV 2 . 





a 


b 


m )it 


fv 


115.74 


-106.73 


23.41 


fo 


10.43 


-12.85 


20.66 


u 


5.50 


-5.07 


14.58 


f- 


-5.36 


4.90 


14.03 



Table 7: Parameters appearing in the fit function (ii) for form factors of the 
B d -> D* d (2010)£u at M 2 = 19 GeV 2 , Mf = 5 GeV 2 . 

the subleading Isgur-Wise form factors for B — > D*£z/ are shown in Table 
8. Note that the values in this Table are obtained with 7\ = T 2 = 2 GeV 
correspond to M 2 = 19 GeV 2 and Mf = 5 GeV 2 which are used in Tables 
[2-7]. 

Table 8 shows a good consistency between our results and the prediction 
of [IS] for Ri at zero recoil limit, y — 1.1 and 1.2, but for the other values 
of y, the changes in present work results are little greater. The values for 
i?2 shows an approximate agreement between two predictions, however the 
changes in the value of R% in our work is also a bit more then [T5]. For both 
Ri and R2, our study and [15] predictions have the same behavior, i.e., Ri 
decreases when the value of y is increased and increasing in the value of y 
causes the increasing in the value of R2. From this Table, we also see that the 
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y 




i i 
i . i 


1 9 


1 3 
1 .o 


1 zL 


1 ^ 

1 . u 


a 2 (RpV 2 \ 

y i \jrc- v j 


1 fiQ 


8 ^7 


fi 4^ 


4 


9 2D 


n D8 


Ri 


1.34 


1.31 


1.25 


1.19 


1.10 


0.95 


R 2 


0.80 


0.99 


1.10 


1.22 


1.30 


1.41 


R 3 


-0.80 


-0.79 


-0.80 


-0.81 


-0.80 


-0.80 




0.50 


0.64 


0.77 


0.94 


1.20 


1.46 


^5 


-0.50 


-0.51 


-0.56 


-0.62 


-0.71 


-0.89 


Ra 


-0.80 


-0.67 


-0.64 


-0.61 


-0.55 


-0.53 


Ri [15J 


1.31 


1.30 


1.29 


1.28 


1.27 


1.26 


R 2 [15J 


0.90 


0.90 


0.91 


0.92 


0.92 


0.93 



Table 8: The values for the Ri and comparison of values with the pre- 
dictions of [T5l. 



i?4 is sensitive to the changes in the value of y. However, the results of R%, 
i?5 and Rq vary slowly with respect to y. Our numerical analysis for l/m^ 
corrections of form factors in Eq. (1251) shows that this correction increase the 
HQET limit of the form factors fy and /+ about 7.1°/o and 6°/o, respectively, 
however it doesn't change the fo and decrease the /_ about 6.5°/o. 

The next step is to calculate the differential decay width in terms of the 
form factors. After some calculations for differential decay rate 



dT 

dq 2 



4™ 2 



7 I Gl I V cb | 2 {{2A, + A 2 q 2 )[\ f v | 2 (4m 2 | 7 | 2 )+ I fo Hi 



Bo 



+ 



167r 4 m 



+ m 



B q 

~r! I 2 



7? 12 



Ba 



m 



111 



<f)) + I fo I 2 
(2m 2 + 2m 2 - q 2 



2 I ' 1 2 

j/ 1 2 m B q I P I /0 _ 2 , o~2 ..2, 



, ,2 m B 9 I P I' 2 



A', 



2 I ' |2 

2 ™ B « 2 (^(/o/+ + /o/- + « - ™\)f + f-)) 
m D * q q 
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2A 



2 I ~? 

m% | p 



2 r 



m 



/ji 2 +K-<fi/;i 2 +? 4 if 12 



+ 2(m 2 - m*)Re{f ft) + 2q 2 f'J* + VK, - m 2 )Re(f + f_ 



(30) 



is obtained, where 



p | = — 



2m^ 

'2 „,2\2; 



A 2 = ^(q 2 -m 2 )(q 2 + 2m 2 )I 
2 g 2 

/o = /o(^D*+m Bg ), 

/v 



fv 



(mn* +m Bg ) 



/L = t § r- (31) 

{m D * +m Bq ) 

The following part presents evaluation of the value of the branching ratio 
of these decays. Taking into account the q 2 dependence of the form factors 
and performing integration over q 2 in the interval mf < q 2 < (rriB q — fn^,*) 2 
and using the total life-times t Bu = 1.638 x 10~ 12 s , r Bd = 1.53 x 10~ 12 s [H] 
and tb s = 1.46 x 10~ 12 s [31], the branching ratios which are the same for 
both fit functions are obtained as: 



B{B S -> D*{2ll2)iv) = (1.89-6.61) x 10 



-2 
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B(B d -> D*(2010)£v) 
B(B U -> D:(2007)£z/) 



(4.36 - 8.94) x 10" 2 , 
(4.57-9.12) x 10~ 2 . 



(32) 



The ranges appearing in the above equations are related to the different 
lepton masses (m e , m M , m T ) as well as the errors in the value of input param- 
eters. Finally, we would like to compare our results of the branching ratios 
with the predictions of CQM model [13] and existing experimental data in 
Table 9. From this Table, we see a good agreement among the phenomeno- 
logical models and the experiment for u and d cases. However for s case our 
results are about 1.7 times smaller than that of the CQM model. Also, there 
is a same behavior between present work results and the experiment. In the 
experiment, the value for branching ratios decreases from u to d. In our 
results also, this value decreases from u to s cases. The order of the branch- 
ing fraction in present work for B s — ► D*iv decay shows that this transition 
could also be detected at LHC in the near future. For the present and fu- 
ture experiments about the semileptonic b — > civ based decays see [32] [37] 
and references therein. The comparison of results from the experiments and 
phenomenological models like QCD sum rules could give useful information 
about the strong interaction inside the D* and its structure. 

In conclusion, the form factors related to the B q — > D*lv decays were 
calculated using QCD sum rules approach. The HQET limit of the form 
factors as well as l/m^ corrections to those limits were also obtained. A 
comparison of the results of form factors in HQET limit with the application 
of the subleading Isgur-Wise form factors at zero recoil limit and others 
values of y was presented. Taking into account the q 2 dependencies of the 
form factors, the total decay width and branching ratio for these decays were 
evaluated. Our results are in good agreement with that of the CQM model 
and existing experimental data. The result of B s — > D*iv case shows a 
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possibility to detect this decay channel at LHC in the near future. 





B s -> D* B iv 


B d -> D* d lv 


B u -> D*Jv 


Present study 


(1.89-6.61) x IO" 2 


(4.36 -8.94) x IO" 2 


(4.57-9.12) x IO" 2 


CQM model 


(7.49 - 7.66) x 10" 2 


(5.9-7.6) x IO- 2 


(5.9-7.6) x IO" 2 


Experiment 




(5.35 ±0.20) x 10~ 2 


(6.5 ±0.5) x IO" 2 



Table 9: Comparison of the branching ratio of the B q — > D*lis 
decays in present study, the CQM model [13] and the experiment [14] . 
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Figure 3: The dependence of f on q 2 at M 2 = 19 GeV 2 , M 2 = 5 GeV 2 , 
s = 35 GeV 2 and s' = 6 GeV 2 . 
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Figure 6: Comparison of fit functions (i) and (ii) for form factor f v for q = s. 
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